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Abstract 

The gAN-y coupling for the A — * N*y decay is calculated in framework of the light 
cone QCD sum rules and is found to be gAN-y = (1>6 ± 0.2) GeV . Using this value 
of 9AN-y we estimate the branching ratio of the A + — ► N'j decay, which is in a very 
good agreement with the experimental result. 
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1 Introduction 



The extraction of the fundamental parameters of hadrons from experimental results requires 
some information about the physics at large distance. Unfortunately such an information 
can not be achieved from the first principles of a fundamental theory. Indeed QCD, which 
is believed to be the candidate of such an underlying theory of the strong interactions, has 
very complicated infrared behavior which makes it impossible to calculate the properties of 
of the hadrons starting from a fundamental QCD Lagrangian. Therefore for determination 
of the parameters of hadrons a reliable nonperturbative approach is needed. Among all 
nonperturbative approaches, QCD sum rules method which was originally proposed by 
Shifman, Vainshtein and Zakharov [|l|] and adopted or extended in many works [[|-||, 
is particularly a powerful one in studying the properties of the low-lying hadrons. In 
traditional QCD sum rules method |l[ the nonperturbative approach is taken into account 
through various condensates in the nontrivial QCD vacuum. 

In this work we employ an alternative approach to the traditional sum rules, namely 
light cone QCD sum rules method, to study A — > decay coupling constant. Light 
cone sum rules is based on the operator product expansion on the light cone, which is an 
expansion over the twists of the operators rather than dimensions in the traditional sum 
rules. The main contribution comes from the lowest twist operator. The matrix elements 
of the nonlocal operators sandwiched between a hadronic state and the vacuum defines the 
hadronic wane function (more about application of light cone QCD sum rules can be found 
in 0-|16] and references therein). 

In general A — > Nj decay is described by the electric quadrapole E2 and magnetic dipole 
Ml transition amplitudes. However, it is well known that the electric quadrapole amplitude 



is very small compared to that of the magnetic dipole amplitude (see and references 
therein). Therefore in this work we consider only the magnetic dipole contribution. 

The coupling constant <7ajv7 is involved in phenomenological models in investigation of 
the many reactions of the strong and electromagnetic interactions and it is expected to 
be measured more precisely in the pion photo production experiments at TJNAL (former 
CEBAF). 

The paper is organized as follows: In section 2 the light cone QCD sum rules for 
the radiative A — > is presented. Section 3 is devoted to the numerical analysis and 
discussion of the results. 



2 The light cone QCD sum rules for A — > ^7 decay 
constant 

In studying the A — > N'y decay constant we first introduce the interpolating currents for 
the A and N baryons H 
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where u, d are up and down quark fields, respectively, and C is the charge conjugation 
operator, a, b and c are the color indices. Note that this choice of nucleon interpolating 
currents is not unique and other choices can be used (see for example ||). 

The overlap amplitudes of the interpolating currents with the baryons are defined as 
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where u A is the Rarita-Schwinger spinor. 

The coupling constant gAN^ for the A — > decay is defined as follows 
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where e p and g p are the polarization vector and momentum of the photon, respectively, e 
is the electric charge. 

According to the QCD sum rules ideology, the quantitative estimates of the (Jan^ cou- 
pling constant can be obtained by equating the two different representations of a suitable 
correlator, written in terms of hadrons and quark-gluon language. For this purpose we 
start our analysis by considering the following correlator 



n(p,g)= /dW^< 7 (g)|T{77 A (0)^(z)}|0> . 
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Saturating (4) by A and nucleons and using Eqs. (2) and (3), we get for the phenomeno- 
logical part of the correlator 
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The main theoretical problem being the calculation of Eq. (4) in QCD. This problem 
can be solved in the deep Euclidean region where p 2 and (p + q) 2 are negative and large. 
After lengthy calculations, at quark level we have obtained the following expression for the 
correlator function 
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4^3 



d A xe lvx 



x • 



"2 {l(q) 
-4( 7 (g) 

-2 (l(q) 



4 (7(g) |d757^47 5 7^7 p 7 5 Tr (<Sl P <S'lp) + ZSlpS'l^^pK 
da a p 

ua a pu 



(6) 



0"«/37p75Tr (S^pS'ifi) + 2S7 p S / 7 At o Q/3 7p75 

2<S7 p 7 5 Tr (757 V 7 P S'7 M ) - 2<S7 p 7 5 7^7m< 5 7p75 + 2^ ip -f p S'^ IJl S^ p ^ 5 



257 P 7 5 Tr ((x Q ^7 pl S / 7 / ,) + 2S f y p <r a pn ll S'yp'y h + 2o Q/3 7 p <S / 7 p <S7 p 7 5 







where S' = CSC = —CSC 1 and iS(x) is the full light quark propagator with both pertur- 
bative and nonperturbative contributions 
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It follows from Eq. (6) that, in order to calculate the correlator function II in QCD, the 
matrix elements (7(g) 1 97a 75 9 1 0) and (7(g) \qcr a pq\ 0) are needed. These matrix elements 
are defined in terms of the photon wave functions as follows flS]— f2~0 
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where the parameter \ is the magnetic susceptibility of the quark condensate and e q is the 
quark charge. In further analysis the path ordered gauge factor 

Vexp (^J dux' 1 A^(ux) 

is omitted since in the fixed point gauge x^A^ = 0. The functions <j>{u) and ip{u) in 
Eq. (8) are the leading twist photon functions, while g\{u) and g7,{u) are the twist-4 
functions. Using Eqs. (6), (7) and (8), and performing Fourier transform for the structure 
[ £ ij,(qp) ~ ( £ p)<Z/Jj we get the following result 
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where P = p + qu. 

As has been noted already, the QCD sum rule is obtained as usual by equating the 
hadronic representation of the correlator (4) with the result of the QCD calculation. In 
order to take into account the contributions of the higher states we invoke the quark-hadron 
duality prescription, i.e., above certain thresholds in and s 2 , the double spectral density 
p(si, s 2 ) for the higher states and continuum coincides with the spectral density calculated 
in QCD. 

After performing double Borel transformation with respect to the variables p 2 and (p + 
q) 2 in Eqs. (4) and (9) to suppress the higher states, we finally get the following sum rules 
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for the AN^y coupling constant 
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where the function 
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is the factor used to subtract the continuum, sq is the continuum threshold and 
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where M 2 and M| are the Borel parameters. Since masses of the proton, and A + are very 
close to each other, we can choose M 2 and Mf to be equal to each other, i.e., M 2 = Mf = 
2M 2 , from which it follows that Uq = 1/2. 



3 Numerical results 



It follows from eq. (10) that the main input parameters of the sum rules are photon wave 
functions. It was shown in |7|, |J that the leading photon wave functions receives only 
small corrections from the higher conformal spin, so they do not deviate much from the 



asymptotic form. Following |18fl-|20fl, we shall use for the photon wave function 
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where u — l—u. The values of the other input parameters we have used: / = 0.028 GeV 2 
and X = -4-4 GeV 2 at the scale (i = 1 GeV, (g 2 G 2 ) = 0.474 GeV 2 and for the 
continuum threshold s we have chosen two different values, i.e., s = 2.8 GeV 2 and s = 
3 GeV 2 . Having fixed the input parameters, one must find the range of values of M 2 for 
which the sum rule (10) is reliable. The lowest value of M 2 is usually determined by the 
condition that the terms proportional to the highest inverse power of the Borel parameter 
stay reasonably small. The upper limit is determined by demanding that the continuum 
and and higher state contribution does not get too large, say less than 30% of the leading 
twist contributions. Both conditions are satisfied in the interval 1 GeV 2 < M 2 < 1.5 GeV 2 . 
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The dependence of the right side of Eq. (10) on M 2 is shown in Fig. (1). It follows from 
this figure in that the above-mentioned working region of M 2 the sum rule is quite stable. 
From this figure one can directly predict that 

^atv 7 AaAtv = (0.0020 ± 0.0004) GeV 5 . (11) 

Dividing this product of couplings by the residues of proton and A currents A at and A a, 
that were calculated in the analysis of mass sum rules for baryons (see also H and |4]]) 

g ANl = (1.6 ± 0.2) GeV' 1 . (12) 

This prediction of the coupling constant permits us to estimate the width of the A — > N'j 
decay. Using the matrix element for the A — > N'-f transition (see eq. (3)), we get for the 
decay width 



2 ( 2 
'9 AN*/ \ m A 



4 m |" v r " A mN 



1 + -^2- (m A - m N ) 2 
3m A 



(13) 



where a is the fine structure constant, m A and are the masses of A and nucleons, 
respectively. Using the predicted value of gAN-y in Eq. (12), the result we get for the decay 
width is 

T ~ 0.65(1 ±0.25) MeV , 

and for the branching ratio of this channel we have (for the total decay width we have used 
Ttot = H3 MeV H) 

r 

B(A -> iV 7 ) = — = 0.0058 . 

±tot 

This prediction is in a very good agreement with the experimental results, i.e., 0.52% < 
B < 0.60% 



In summary, we have calculated AA^7 coupling using the light cone QCD sum rules. 
Our prediction on the branching ratio is in a good agreement with the experimental results. 
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Figure captions 

Fig. 1 The dependence of gANy^A^N on the Borel parameter M 2 . In this figure the solid 
and dash-dotted lines correspond to the threshold values s = 2.8 GeV 2 and s = 3.0 GeV 2 , 
respectively. 
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